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CN Abstract 



We discuss exotic properties of charged hydrodynamical systems, in the broken 
superconducting phase, probed by electromagnetic waves. Motivated by general 
arguments from hydrodynamics, we observe that negative refraction, namely the 
propagation in opposite directions of the phase velocities and of the energy flux, is 
expected for low enough frequencies. We corroborate this general idea by analyzing 
^ a holographic superconductor in the AdS/CFT correspondence, where the response 

functions can be explicitly computed. We study the dual gravitational theory both 
in the probe and in the backreacted case. We find that, while in the first case the 
refractive index is positive at every frequency, in the second case there is negative 
refraction at low enough frequencies. This is in agreement with hydrodynamic 
considerations. 
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Introduction 

In the last years there have been much progress in realizing artificial materials which 
exhibit exotic electromagnetic phenomena. These artificial materials, commonly called 
electromagnetic metamaterials, consist of periodic arrays of metallic rods Ql and split ring 
resonators (2). Among others, one of the most fascinating properties is the possibility of 
negative refraction, i.e. the energy flux of the electromagnetic wave flows in the opposite 
direction with respect to the phase velocity. This has been theoretically suggested by 
Veselago in 1968 [3] and only recently experimentally realized (4j|5). The main difficulties 
in the realization of such systems is the necessity of strong magnetic response, and their 
high degree of electromagnetic dissipation. 

In (6) it was argued that negative refraction is a common feature for transverse elec- 
tromagnetic waves probing relativistic charged hydrodynamical systems at low enough 
frequencies. This idea was supported by an explicit computation of the permittivity and 
the permeability of a strongly coupled system that admits a dual gravitational description. 
This system is supposed to describe a four dimensional charged strongly coupled plasma. 
Moreover in [7] the same topic was investigated for a charged black hole in four dimen- 
sions, and in [8j[9] for systems with D7 flavor branes. Especially this last development 
seems promising for more realistic applications to the Quark Gluon Plasma physics. 

One may wonder about the electromagnetic response properties at low frequencies 
of other materials IIOJ. For example in the last years an interesting proposal of super- 
conductors as NIR (negative index of refraction) materials or metamaterials have been 



investigated (for review see 11 and references therein). It was shown that in the radio, 
microwave and low-terahertz frequency range the superconductors can behave as meta- 
materials and exhibit negative refraction. The interest in building NIR superconductors 
is triggered by the fact that they can reduce the dissipation (lossy) of the electromagnetic 
energy, thus evading the constraints on the performance of the usual metamaterials and 
providing a convenient setup for experimental investigations. 

In this paper, motivated by the experimental developments, we investigate the refrac- 
tive index of a class of superconductors in which the response functions are rather simple 
and their dependence on the transport coefficients is known. They are superconductors 
that at low frequencies and wave- vectors can be described by hydrodynamics. We observe 
that quite generically it is possible to find frequency ranges where these media behave as 
NIR materials, when opportune choices of the transport coefficients are made. 

Then we verify this general hydrodynamic prediction in a class of superconductors in 
which the transport coefficients can be explicitly computed. The class of materials that 
we are describing are referred in the recent literature as holographic superconductors. 
They are associated to strongly coupled gauge theories dual to some weakly coupled 
gravitational system. The latter weakly coupled theory simplifies the computation of 
the correlation functions describing the linear response of the material to an external 
electromagnetic source. 

The possibility of describing a superconductor in the context of the gauge/gravity 



duality was shown for the first time in 12-14 . According to this duality, the conformal 



field theory describing the quantum critical region is mapped to a gravitational system 
in one higher dimension whose asymptotic spacetime is AdS. The U(l) symmetry of the 
field theory is associated to a gauge field living in the bulk" . The phase transition in the 
gravity theory is set by a scalar field which breaks the U(l) symmetry below a critical 
temperature. This is the mechanism first explained in [12]. Strictly speaking, the field 
theory U(l) symmetry is a global one, so we should talk about superfluidity rather than 
superconductivity. Anyway, a weak gauging of the symmetry is consistent, and it opens 
up the way to the description of a superconductor in AdS/CFT. In particular, below the 



critical temperature the breakdown of the U(l) adds a pole in the conductivity 13, 14 



thus modifying the electric permittivity. Many examples of this type have been studied, 



and we refer the reader to the reviews 15 and references therein. 

A previous study of the refractive index in holographic superconductors was done 
in 



16 where the authors observed that the refractive index is positive for every frequency. 
That result holds for a holographic superconductor in a specific limit, called the probe 
limit. This corresponds to the limit in which the metric background does not fluctuate, 
and the response functions are computed only from the fluctuation of the Maxwell field. 
When the background is fixed the correlation function of the transverse electromagnetic 
current, below the critical temperature, is dominated by the Goldstone modejj Then the 
conductivity only has the pole due to the Goldstone boson which is not enough, at least 
at the leading order for low frequencies and wave-vectors, to obtain a negative refractive 



3 This result can be also understood quite generally from the hydrodynamical equations. 



index. 

Allowing background fluctuations corresponds in the hydrodynamical description to 
also considering the stress-energy tensor conservation equations, which imply translational 
invariance in the field theory due to momentum conservation. In this case there is also a 
diffusive pole in the current correlator. At low frequencies this suggests the existence of 
negative refract ion jj 

In this paper we look for negative refractive index in backreacted holographic super- 
conductors. We indeed find that if the backreaction is taken into account the refractive 
index becomes negative at low frequencies. This result corroborates the general claim that 
in systems described by hydrodynamics at finite temperature and finite charge density the 
refractive index becomes negative in the low frequencies regime, because the backreaction 
is proportional to the charge density. 

The rest of the paper is organized as follows. In section [T] we review the computations 
of the linear response of media and the electromagnetic response properties, focusing on 
the negative refractive index. Then in section [2] we illustrate the general prediction on the 
refractive index of media described by hydrodynamics at finite temperature and charge 
density. In particular, we observe that negative refraction is expected for a superconductor 
if the system has a finite charge density and the diffusive contribution cannot be neglected 
with respect to the Goldstone mode typical of the superconducting broken phase. In 
section [3] we introduce the bulk setting for the study of a holographic superconductor. In 
section [4] we review the probe limit by showing even at the analytical level that in that 
case the refractive index is always positive. In section [5] we allow the backreaction of the 
matter fields on the metric. In this case we observe that increasing the backreaction the 
refractive index becomes negative at low frequencies. Finally we conclude. 

1 Linear response functions and negative refraction 

In this section we review the formalism of the linear response theory of a continuous 
medium to the EM field, and its application to the study of the negative refractive index. 
The electrodynamics response of a medium is described by the electric and magnetic 
fields and inductions E,B,D,H. These fields are related by the constitutive equations (in 
Fourier space) 

Diiw) = e ij (w)E j (w) , Bi(w) = Hij^H^w) (1) 

The response functions e and /i are the permittivity and permeability of the medium. They 
are complex quantities and depend generically on the frequency of the electromagnetic 
wave. 

In hydrodynamical systems, which the strong coupling media we are going to describe 
are a particular example of, non local effects can be relevant, and hence the response 
functions are in this case functions also of the wave vector k. This phenomenon is known 



6 We expect that breaking the translational invariance with small impurities does not qualitatively 
alter this conclusion. However a more detailed study would be interesting. 



as spatial dispersion. The magnetic and electric response of the system can be here just 
described in terms of the field D, B and E (171, where 



D t (w, k) = €ij(w, k)Ej(w, k) (2) 

If the medium is isotropic the permittivity e^ is decomposed in a transverse and in a 
longitudinal part, ex and ex,, with 

€ij{w, k) = -r^tLiw, k) + I bij p- I e T (w, k) (3) 

The EM properties of the medium are analyzed by solving the dispersion relation, obtained 
by the Maxwell equations. Here we focus on the transverse propagation, which has the 
dispersion relation 

k 2 
e T (w,k) = — = n 2 (4) 

w z 

where n is the refractive index of the medium. It is a complex quantity which, once we 

impose the dispersion relation, is a function only of the frequency. The real part of n is 

the usual refractive index, while the complex part encodes the dissipation. 



If the spatial dispersion is small, we can expand the transverse permittivity as 18 



k 2 ( 1 \ 

e T (w, k) = e(w) + _ l ___+.. . (5) 

ur \ H(w) J 

This expansion connects the EBD approach with the EBDH onej^j Indeed the expression 
(Jsl) has been arranged in such a way that if we impose the dispersion relation Q we obtain 
the usual dispersion relation for an electromagnetic wave in the e - fi approach, that is 
e/i = n 2 . 



It has been shown 19 that in case of dissipative media described by complex e and jjt, 



the sign of the refractive index coincide with the sign of the function 

n DL {w) = Re(e(w))\n(w) | + Re(ji{w))\e(w) | (6) 

This index was originally derived for dissipative systems without spatial dispersion. It 
corresponds to the requirement that the orientation of the Poynting vector S is opposite to 
the orientation of the phase velocity which is Sign[Re[n]]. In the case of the propagation 
of a transverse wave the Poynting vector can be written as 

S = *g)W (T) 

Imposing Sgn(i?e(n//i)) =-Sgn(i?e(n)) is equivalent to require that ([6]) is negative. Here, 
in the presence of spatial dispersion the definition of the Poynting vector is slightly mod- 
ified. Anyway, as shown in [6] , the relation (IT]) is still valid if the spatial dispersion is not 

7 However it is important to observe that the function fi(w) we introduce here is an effective magnetic 
permeability, that contains the magnetic response of the medium 17 plus some k 2 corrections to the 



electric response. In the case of dispersive media it is the right quantity to study waves propagation 18 . 



too large and uol is still a measure of the sign of the refractive index. In the rest of the 
paper we will use ^ as an index to study at which frequencies the system at hand shows 
negative refraction. 

The response functions e(w) and fi(w) can be computed from the linear response 
theory applying an external field Aj. Indeed the resulting U(l) current Jj is proportional 
to Aj through the relation Jj = q 2 GijAj, where q is the electric charge and Gij is the 
retarded correlator of U(l) currentr] For an isotropic medium, we decompose G in a 
similar fashion to pL The transverse component Gt(w, k) is related to e^ byn 



47T 



e T (w, k) = H T^q Gt(w, k) 



As anticipated above in the case of small spatial dispersion we can expand this expression 
at the second order in the wave- vector (|5). Denoting Gt(w, k) = G T (w) + k 2 G T \w) + 
0(k 4 ), the permittivity and permeability read 



Air 



e{ w ) = 1 + —q 2 G^(w 



2/?(0), 



w 



^ = A o^(2), \ ( 9 ) 

1 -Aiiq 2 Gy(w) 

In particular the conductivity is defined as: a(w) = —iq 2 G T \w)/w and e(w) = 1 + 
47ria(w)/w. 

2 Hydrodynamics and negative refraction 

Before doing any explicit computations it is worth noting that we can infer the generic 
behavior of the refractive index for the type of systems we are going to consider in the 
rest of the paper from hydrodynamic arguments. As explained in the previous section, 
to reach this goal we need to know the generic form of the retarded correlator of the 
transverse current. This goal is actually possible at least for the leading order in the low 
frequencies w and low wave vector k expansion. Indeed in the limit of low frequencies 
and long wave lengths a system can be typically described by hydrodynamic equations. 
These equations describe the long time dynamics of the effective macroscopic degrees of 
freedom of the system: conserved charge densities and phases of order parameters. The 
linearized hydrodynamic equations describe the response of the macroscopic degrees of 
freedom to the system due to the application of a small external field. 

In this section we will quickly review the general argument in the specific case of the 
relativistic superconductor, the interested reader is referred to [To] for more details. To 
obtain the correlator of the transverse current J T it is useful to study the dynamics of 



8 In many cases the explicit q 2 dependence in this relation is absorbed in the Green function. Here we 
prefer keeping it explicitly to stress that the EM field in AdS/CFT is introduced as an external field. 

9 Here we follow the opposite convention of |6| f or the sign of the Green function to be consistent with 
the conventions common in the literature (es. [20 ). This changes the sign of the q 2 contribution in ct- 



the transverse momentum density tt . As usual in the hydro dynamic regime 20,21 the 
transverse part of the momentum density decouples from the other degrees of freedom 
and it satisfies the diffusive equation: 

(0 t - T?V 2 )vr T = (10) 

where V = (e+P ^_ s, and r\ is the shear viscosity, e the energy density, P the pressure, 

p s the charge density of the superconducting part of the fluid and p p the electrostatic 
chemical potential. 

On the other hand linear response theory gives the non equilibrium momentum density 
of the system as a function of its retarded correlator computed in the equilibrium state 
as: 

iw G^-t^O, k) 

where n T (t = 0, k) is the Fourier transform of the initial value for the momentum density. 
Linear hydrodynamics and linear response theory at long time and large scale are 



two different description for the same physics, and they must coincide 20 21 . If we 



equate (11) with the solution of the initial value problem (10), taking the k — >■ limit for 
the correlation function and considering that G„-t(0, 0) = (e + P — p s p p ), we obtain the 

classical expression for the momentum correlator: G 7T t(w, k) = _ jw lp fc 2 ■ 

We can now compute the retarded correlator of the transverse current: Gjt(w, k). The 
induced transverse current J T is the response of the system to an external electromagnetic 
field v4g Xi . In presence of an electromagnetic field the transverse velocity field of the system 
is no more proportional to tt t but to ir T — A T p t , due to the minimal coupling of the system 
to the vector potential, where p t = p n + p s is the total charge density: the sum of the 
charge density of the normal fluid plus the charge density of the superfluid component. 



The transverse current is proportional to the transverse field velocity of the system 20 - 23 
and in particular: 

J T = PnV T = T ~^ S - T ~ — |^ S A T (12) 

(e + P - psPp) (e + P - p s p p ) 

We then obtain Gjt(w, k), if we assume that all the electromagnetic fields in the system 
are external: namely A = A ext . This hypothesis can be a bit strong and amounts to 
treating even the internally generated electromagnetic field as an external perturbation. 
We believe that it can be justified in the case we are studying in the main text, where 
the system itself is dominated by the non-abelian strong interactions and there is no 
dynamical electromagnetic field if we do not couple it with an external field. 



Using equation (12), the fact that J T = GjrA^ xt and the explicit form for G^t^w, k), 
we obtain: 

GjT f w fc N = %W {e+P-Psto) Ps(pt-Ps) _ JWB c ^ 

-iw+ {e+P l Pallp) k 2 (e + P-p s p p ) -iw + Vk 2 



where £>, C and T> are real constants. In the limit of normal fluid p s — > 0, p n — > pt, the 



constant contribution C to Gjt goes to zero, and (13) reduces to the correlation function 
discussed in (6|, for a normal relativistic, translation invariant, isotropic system at finite 
charge density. 



From the generic form of (13) we can obtain G(°\w) and G^ 2 \w) and hence the generic 
form of e(w) and p{w): 

e(w) = l-^(B + C), p(w)= i _ 4 1 (14) 

^ in 



«- 



and of the conductivity: 

a( w ) = l g 2 (£ + C) (15) 

w 

from which it becomes clear that the C is responsible for the infinite DC conductivity of 
the superconducting phase, while B is responsible of the infinite DC conductivity of the 
normal phase, associated to the translation invariance of the system. 



From the equations (14) one can verify that for low enough frequencies udl(w) is 
negative. 

Hence we conclude that an isotropic and homogeneous superconductor at equilibrium, 
with relativistic invariance and at finite charge density, has a negative refractive index for 



10 



low enough frequencies in the hydrodynamical regime. 

Differently from system in the normal phase, in the case of superconductors the imag- 
inary pole in the conductivity has two different contributions; one is associated to the 
translational invariance and it comes from the charge density of the normal fluid p n , and 
the other is associated to the superconducting phase of the system and it comes from p s . 
Both terms contribute to make Ke(e(w)) negative at low frequencies. This fact is usually 
a necessary condition to obtain negative refraction. Anyway also the imaginary part of 
the permeability conspires to give negative refraction. This suggests that the diffusive 
pole is crucial in order to have a negative refractive index. Indeed, by varying the amount 
of backreaction, we will be able to explore the interplay between the two poles and to 
show that by raising the size of the backreaction the refractive index becomes negative at 
low frequency and within our approximation. 

Moreover in (6J, even if negative refraction was found, the system was shown to be 
highly dissipative. In superconductors the real part of the conductivity has a gap for 
low enough frequencies in which the conductivity is almost zero, this implies that the 
imaginary part of the electric permittivity is very small and hence the dissipative effects 
are smaller than for a normal fluid. Indeed, due to the presence of a physical pole in the 
imaginary part of the conductivity and a gap in the real part of the conductivity, the 
superconductors can be convenient systems with negative refraction and low dissipation, 
and they have been already used to obtain negative refractive metamaterials in (ill . 



'Actually the conclusion is more general and can be extended to system which are not relativistic 



invariant and are not in the superconducting phase. See 10 . 



3 Holographic superconductors: the bulk setting 

In this section we start the study of the basic setup for the analysis of the superconduc- 
tivity of a medium through the techniques of the gauge gravity duality The minimal 
ingredients to build a holographic theory of s-wave superconductivity are a bulk abelian 
gauge field A^ and a charged massive scalar field <fi coupled to gravity. We are then led 



to consider the following action 13[ 14 

1 

s ■ • — 



d £-*/—# 



^yR+^-l^F^-Wf-m 2 ^ 



(16) 



where the Newton constant G5 is related to the backreaction by k 2 = 87rG 5 



fit/ 



d»A v 



D, 



-d v A^ 



(17) 



are the gauge field strength and the covariant derivative respectively, V M is the metric 
covariant derivative and q s is the charge of the scalar field in five dimensions. Note that 
we used the freedom of rescaling the gauge field so that its kinetic term is fixed, and q s is 
the charge of the scalar field in these units. By also defining the real current 



J^ = i 



V DA - D„ 



:u 



the Einstein, Maxwell and scalar equations of motion resulting from (16) can be written 
as 

6 



G 



i-"' 



L 2 



9^v 



k z T, 



\o/ 



V^F, 



i-"' 



(D^D^ - m 2 ) 







(19) 

(20) 
(21) 



with G^y the Einstein tensor and T^ the stress-energy tensor. 



Throughout this paper we only consider the plane-symmetric ansatz for the metric 



ds z 



-r 2 f{r)e 2v{r) dt 2 + 



dr 2 



r 2 f(r] 



+ r 2 dy^ 



(22) 



where the condition /(r^) = determines the radius r^ of the black hole. Moreover, 
the equation of motion for / sets f(r) — > 1, i.e. the spacetime is asymptotically AdS: 
according to the AdS/CFT, the conformal symmetry is recovered in the UV limit of the 
dual field theory. 

We take the following ansatz for the gauge and scalar fields 



At 



h{r) 
Mr) 



A = 



r,x,y,z 



(23) 



which corresponds to a homogeneous system. Because the scalar field is supposed to be 



charged under the bulk U(l) gauge symmetry, we wrote it as a complex field in (16). 
However, the r component of the Maxwell equations sets its phase to a constant, and 
because this constant does not affect the other equations, we set it zero without any loss 
of generality. Note that the phase of the scalar field cannot be always set to zero. 

r 2 

We find more convenient to change the radial coordinate to u = -| so that it is a 
compact direction: the horizon is located at u — 1 and the boundary at u — 0. In this 



coordinate system, the ansatz (22) becomes 

f{u)eMU) dt 2 I rl ^ I - 2 
u 4u 2 f(u) u 



ds 2 = - J -^^ dt 2 + ^-^- + '-±dy? (24) 



Unless otherwise specified, from now on we use the coordinate u and we set r^ = 1. 
The Hawking temperature of the black hole is given by 

_ /'(l)e»W 

Th ~ ~^r~ (25) 

and is identified with the dual field theory temperature only if i/(0) = 0. 

4 The probe limit 

In this section we review the refractive index of the holographic theory described above 
in the limit of k = 0. Even if this analysis was already performed in n6\ here we prefer to 
start from this simplified example to fix the notations which will become important for 
the backreacted case. 

We start with the analysis of the equations of motion and explain the derivation of 
the Green functions in this limit. Then we study the numerical solutions and we observe 
that the system has a positive refractive index even at low frequencies. We conclude this 
section comparing the numerical results with the analytic computation. 

4.1 The holographic setup 

In the probe limit k 2 = the Einstein equations decouple from the gauge-scalar sector 
and the analytic solution for the metric is the usual Schwarschild-AdS black hole: The v 
function is constrained to vanish, while f(u) = 1 — u 2 . We are left with a system of two 
ordinary differential equations for the background 

Because they are non-linear, an analytic solution where both fields are nontrivial is very 
hard to find. However, mainly motivated by the corresponding numerical solution, an 

9 



analytic perturbative solution can be computed by assuming a small scalar field 24 . We 



will use this method later, when we will compare our numerical computation with the 
perturbative analytic solution to check the strength of our numerical algorithm. 



The system of equations (26) has to be supplemented by boundary conditions for the 
solutions to be completely determined. At the horizon u = 1 we require the gauge field 



to have a finite norm and the equations (26) to be regular there. This amounts to set 



h{l) = 



4>'iX) 



m 



/'(I) 



;i) 



At the boundary, the general asymptotic behavior of the fields is 



h(u) = fip- p t u + 



u 



(D u A/2 + 



(2) M 2-A/2 



+ ... 



(27) 



(2* 



where \i v and p t are the field theory chemical potential and total charge density, respec 



25 



tively, and A == 2 ± vi + m 2 . As long as A > 1 we can choose both signs for A 

for concreteness, we will choose only the upper one in what follows. Thus, we have a 

normalizable solution only when 

(2) = (29) 

According to the AdS/CFT dictionary, cf)^ 1 is interpreted as the expectation value (O) of a 
dual field theory operator with conformal dimension A, and <ffi-> is the value of its source. 



The boundary condition (29) corresponds to the spontaneous breaking of a symmetry in 



the dual field theory: the charged operator O is acquiring a vacuum expectation value 



read from the general falloff (28). 



without being sourced by any field. Once (29) has been imposed, the other quantities are 



In the probe limit the Hawking temperature is a constant 



H 



1 

An 



and an analytic solution of (26) is 



h[u) = p p {l — u) 



(30) 



(31) 



with a vanishing scalar field. The equation (31) represents the normal phase of the dual 



field theory, i.e. the non-superconducting phase above the critical temperature. 

The superconducting phase corresponds to a solution in which the scalar field assumes 
a nontrivial profile. This scalar field is identified with the condensate triggering the phase 
transition. By fixing a reference scale with the chemical potential /i p one can look at 
the evolution of y^)//i p as a function of T / p p . Once the ratio T/p p is lower than a 
certain critical value, the normal phase solution (31) produces an instability towards the 



formation of a hairy black hole. Indeed, one can see that the scalar effective mass in (26) 



is proportional to 12 



m e// = m 2 - 



Qs h ( u ) 



-u 



m 2 + q 2 s g u (u)h(u) 2 



(32) 



10 



When the negative second term is such that m 2 ^ falls below the BF bound |26|, fn 2 ^ < 
m 2 BF = —4, for a range of the coordinate u, the instability occurs and the solution with 
the vanishing scalar field is no longer the vacuum of the theory. 

Till now, our discussion has been quite general. From now on we specialize to the case 
m 2 = -3 with A = 3. 

The next step is to compute the retarded correlator for the transverse current in the 
superconducting phase. To reach this goal we need to study the oscillations of the five 
dimensional gauge field SA x (u, t, z) = a x (u)e~ twt+lkz transverse to k, in the superconduct- 
ing five dimensional gravitational background. As in the standard approach we define 
the dimensionless frequency and wave-vector by normalizing w and k by the temperature 
and then we work with the dimensionless quantities w and k. Due to the isotropy of the 
system we chose the x component of the field, but the same is true for the y component, 
they are indeed both transverse to k, that in this particular case is along the z direction. 

Following the general holographic dictionary |27| we can then compute the GiJ(w) 

(2) 

and Gxd(w) contributions to the current-current correlation function from holography 

and from that obtain the necessary information on the refractive index in the case of 

small spatial dispersion. 

In the probe limit it suffices to allow fluctuations of the x component of the Maxwell 

field, because all the other equations decouple. The linearized equation for this fluctuation 

is 

f'(u) , 1 / w 2 k 2 q^(uY 



a"( u ) + ^(u)---^ + —-- - — - a x (u) = (33) 

}{u) 2uf{u) \2f(u) 2 u 



■tw 



Near the horizon we impose the fluctuation a x to be proportional to (1 — u)~^ , where the 
minus sign corresponds to the ingoing boundary condition. Near the boundary, u — > 0, 
the solution behaves as 

a x (u) = a X)0 + a x>1 u - a Xi0 (w 2 - k 2 )u log(w) (34) 

where a X fl and a X) i are integration constants. 

By using Lorentz invariance, and because we only need up to the k 2 coefficient of 
the Green function, we expand a x (u, w,k) = a x (u,w) + k 2 ai (u,w). Then, the near 



boundary solution (34) reads 

4 0) = 4 J + 4°N-w 2 4>iog( M ) (35) 

4 2) = a® + a S u + a x °l u log(u) (36) 

The Green function is then 

a (0) a (0) fa {2) a {2) \ 

Gi°2(w) = ^ + cw 2 , G£>(w) = ^ (-^--^ )-c (37) 

a x,0 a x,0 \ a x,l a x,0/ 

where c is an arbitrary real constant fixed by requiring e(w — > oo) — > 1, i.e. the system 
at large frequencies behaves like the vacuum. 

11 



In 16 this system was numerically solved and it was observed that nm > in the 
whole frequency range. In the next section we will confirm this result when the effect of 
the backreaction are negligible. In the next subsection, we show that near T c this result 



can be recovered by looking at the analytical computation of 24 



4.2 Connection with the analytical results 

The numerical results that we discussed in the probe limit show that for every temper- 
ature below T c the refractive index is a positive quantity. This is consistent with the 
observation that at T > T c the U(l) symmetry is unbroken and the physics is described 



by 28 . The T ~ T c behavior can be explained even at analytical level. We refer to 
the system studied in [24] as the analytical example to understand the physics of the 
superconductor near T c . The relevant details necessary to compute the Green function of 



an holographic superconductor at non zero wave-vector k can be found in 24 . Here we 



show the correlator G xx at higher orders referring to 24 for the details on the setup. 

The analytical transverse Green function is studied for low values of the frequency w, of 
the wave- vector k and of the vev of the condensate, which we assume to be £ ~ (1 — T/T c ). 
These low values of the condensate are associated to temperatures near T c . 

The refractive index at T ~ T c is obtained by the study of the Green function in this 
regime. By computing some higher order correction with respect to [24] we obtain 

1 *£ 2 w e m 2 wk 2 | g 2 k 2 (7r 2 -6»7rlog2 + log 2 2) 

Gt w, k, t) = «w — w log 2 (38) 

y ' ' w 2 & 16 4 32 64 y J 

The leading contribution to the electric permittivity at low frequencies is computed from 



(38) by applying rt9J) and it is 



<w) = l + W(4i+fl-SU-¥') (39) 



4 w 2 \ 16/ w 

while the magnetic permeability does not have any pole in w and it can be well approxi- 
mated for small q by 

2 

//(w) = 1 + ^- (£ 2 (tt 2 - 6irr log 2 + log 2 2) - in 2 w) (40) 

We conclude that the leading contribution to the index tidl is 

{0 if w < y/Trq£ 

/ 9 rfA (41) 

2M-^_|-J otherwise 

For lower temperatures the approximation of small £ does not hold, and the numerical 
computation is necessary. As observed in the last subsection even at lower temperatures 
the refractive index remains positive. 

12 



5 The backreacted case 

In this section we study the refractive index of the holographic superconductor by taking 
into account the effects of the backreaction. 

In this case the Goldstone mode in the superconducting phase introduces a gap in the 



real part 29 and a pole in the imaginary part of the electric conductivity a, and they 



depend on the Newton constant, as first noted in 30 31 . Moreover this usual pole of 
the superconductivity, related to the U(l) breaking, mixes with a new pole due to the 
translational invariance. Indeed, differently from the probe limit, where the background 
was fixed, and translational invariance was broken, here the fluctuation of the background 
are allowed. This restores the translational invariance and generates a pole in the Green 
function, because an external electric field uniformly accelerates the charges, that cannot 
relax due to momentum conservation. The EM properties of this pole were already studied 
in [6], and at low frequencies negative refraction was shown to be generic. For low k the 
dominating pole is the superconducting one and the refractive index is positive, as in the 
probe limit, while for higher values of the backreaction the negative refraction is allowed 
as predicted from hydrodynamics. 

5.1 Holography 

In the backreacted case the Einstein equations do not decouple from the gauge-scalar 
sector anymore, because k^O. An analytic solution for the metric is a RN AdS black 
hole in the unbroken phase 

,, N , 9N 2K 2 h(l) 2 u 3 ( 1\ 
f(u) = (1 - u 2 ) + ^— (l - -J (42) 

The analysis is completely analogous to that of section El so we will only sketch it here. 
The background system of equations is 

2 3i/(«) qlh{u) 2 <P{u) 2 
HU) + ^^ + 4eMu) uf{u) 2 = ° ( 43 ) 

mV(«) + ^(nmn) _ <m + /WOO + ^^ + 0>) = Q {u) 



4u 2 f(u) 4e 2 »Wuf(u) 2 u f(u) 

h(u)-h {u)v (u) - 4m2/(m) = (45) 

_ m 2 K 2 <P(u) 2 _ q 2 s UK 2 h(u) 2 <P(u) 2 + 2u 3 k 2 h! \uf f (u) 

[ JW) 3 3eM«)/( w ) + l4bj 

+uf{u) _^fWW^ 



where the fields satisfy equations (27)-(29) and /(l) = 0, u(0) = 0. Once a solution is 



found, the field theory temperature is computed by (25). 
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The response of the system is obtained from the correlation function. As usual these 
are computed from the holographic perspective in terms of the boundary values of the 
fluctuations of the Maxwell field and of the metric. We observe that the only fluctuations 
necessary to obtain the transverse current correlation function Gj x j x = G xx are a x and 
the metric components g x t and g xz . Moreover we rise the index x by using the metric 
component gfi*, obtaining gf = g% x g x t = ug xt and g x z = g xx g X z = ug xz to simplify the 
equations. The final set of equations for the fluctuations is 

wg' t x {u) + e 2v{u) k /(«) g' z x (u) + u k 2 w ti{u) a x (u) = (47) 

w(wrf(«) + k£(«)) (I f'(u) Au) \ g > :{u)+gT{u)=Q (48) 



Ae 2u uf(u) \u f(u) 

k(kgf(u)+wg x (u)) qlK 2 a x {u)h{u)4>{uf 2 , , 

+ (l + v'(u)\g' t x (u)-g'r(u) = (49) 



-a x {u)+ ~,,. g t {u) + 



Au 2 f(u) 2 K ' e 2v f{u)- 

+ ( y^ + ' / («))<(«) + <(«)=o ( 5 °) 



where the equations (47), (|48j) and (49) are not independent. Once we derived the equa- 
tions of motion we can follow the procedure explained in the probe limit case to compute 
the Green functions numerically. The difference is that in this case we have a set of cou- 
pled equations and more care has to be taken in the expansion. Moreover in this case the 
relevant contribution to the boundary action is 

S Bd . = / ( ~^a x (u) (f(u)a' x (u) - pg x (u)) + _— #f ( w ) (ug ,x (u) - g x (u)) 

Jin) g z t (u)(ug' z x (u)-g x z (u))Y ' (51) 



2^ 2 

The second term plays a crucial role in the backreacted case because it leads to the 
diffusive pole and to the appearance of the negative refraction in the low frequencies 
regime. 

5.2 Prelude: numerical computation of the Green functions in 
the backreacted case 

According to the prescription for computing the Minkowski space correlators [27|, the 
Green functions are given by the second derivative of the on-shell boundary action with 
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respect to the boundary values of the dual gravity fluctuations. Numerically, this pro- 
cedure is very difficult to implement. Thus, we will extract the Green functions in the 
following way. 

Let us first note that, near the boundary u = e, we have 



sfOO 



9l\e) 



aJ£) 



A x txt 

e 



A, 



(G xtxt #(0) - G xtxz g*M + G xtx a,(0)) 
(-G xtxz gt(0) + G xzxz g x z (0) - G xzx a x (0)) 



(52) 



A, 



((G^x ~ A xx ) g*(0) - G xzx g x z (0) + G xx a x (0)) 



where we discarded all the terms, including the logarithmic ones, which will be subtracted 
by the holographic renormalization procedure, and do not affect the present discussion. 
They will be taken into account by adding a constant to the Green functionsp] The 
coefficients A^ are read from the boundary action (51); for instance, A xtxt = (2uk 2 ^ 



The equations (52) are a consequence of the on-shell boundary action (51). Given the 



above equations, we conclude that, for instance, 

d 2 (gt'(e)gm) 
dgf(0)* 



G 



xtxz 



A, 



(53) 



We now consider the linearized set (47)-(50). Once we have solved the set (48)-(50) 



we still have to satisfy (47). We interpret it as a set of constraints on the different Green 
functions, in the following sense. If we multiply, for instance, the equation by gf(u), take 



the u — ¥ limit, take the second derivative with respect to gf(0) and use (53) and its 
analogous, we arrive at 

w G xtxt + k G xtxz = (54) 



where we used A 



xtxt 



-A 



argument applies when we multiply by g x (u) or a x (u). 



Equation (|54|) is nothing than a Ward identity. A similar 

Thus, we have a total of six 
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equations, namely the set (48)-(50), plus the three constraints above, for six unknown 
Green functions (assuming the Gij = Gji symmetry). Because we are interested in the k° 
and k 2 coefficients of the Green functions, we expand both the functions gf(u),g z (u) 



a T \u\ 



and the correlators in powers of k. Using again (47) and Lorentz invariance, it is easy 
to convince oneself that g x (u) and a x (u) only contain even powers of k, while g x (u) only 
contains odd powers of k. Accordingly, the correlators which do not involve the coordinate 
z, together with G xzxz , only contain even powers of the momentum, while the expansion 
of the other two correlators starts with a linear term in k. The last statement is obvious if 



one looks at equation (54), or reminds the underlying Lorentz invariance. This procedure 



allows us to numerically compute the desired Green functions even when more than a 



11 This is very usual in holographic renormalization. The arbitrary constant corresponds to an ambi- 
guity in the choice of the holographic renormalization scale. We will fix this constant by the physical 
requirement that the real part of the electric permittivity approa ches 1 at large frequencies. 

12 Note that when we multiply by a x (u), also the third term in (47 1 contributes. 
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fluctuation is involved, and we are not able to decouple the equations. To have further 



support of the validity of this procedure we checked it against some analytic results 28 
In the next subsection we will use it to show that our backreacted model exhibits negative 
refraction. Finally, note that, because in the probe limit the entire procedure reduces to 
the usual definition of the single-field prescription, it is perfectly consistent with previous 



results 16 



5.3 Numerical results 

The following step consists in numerically solving the coupled set of equations of motion 



for a x (u), g x (u) and g x (u) and extracting the Minkowski correlator 27 from the boundary 
action by using the procedure explained above. After we obtain the numerical behavior 
for the current correlator G xx (w,\s) we can discuss the EM properties of the backreacted 
holographic superconductor. We start by showing the different results for the permittivity 
and the permeability at different values of k and T/T c . For the backreaction we choose 
k 2 = {1(T 10 , .1; .2; .3} while for the ratio T/T c we choose T/T c = {.45; .6; .75; .9}. We plot 
the real and imaginary part of the electric permittivity and magnetic permeability. 

We observe from the figures [T]j4] that the real part of the permittivity is negative at 
low frequencies in all the cases. This result is expected both in the k — > limit and in 
the fully backreacted case. Indeed in the first case the Goldstone mode dominates over 
the diffusive pole and this results coincides with the probe limit. In the case of a fully 
backreacted superconductor this contribution is corrected by the presence of the diffusive 
pole at zero wave- vector k. The real part of the permeability is always positive and it is 
more peaked near the origin as far as the backreaction is raised. The imaginary part of 
the permittivity is always positive and has a pole at zero frequency. The imaginary part of 
the permeability is negative if the backreaction is very small, q] When the backreaction is 
turned on, the diffusive pole of the green function introduces a pole in the imaginary part 
of G^ and thus in the imaginary part of the permeability. These results are perfectly 



consistent with the general prevision from hydrodynamics ( 14 ) . 

In figure [5] we plot the behavior of the index udl until w = 2, while in figure [6] we 
zoom the plot to the frequency region where negative refraction is sizeable. When the 
backreaction is turned on, a region of negative refraction opens up at low frequencies. 
The size of the region increases with increasing the backreaction parameter. This con- 
firms the fact that the negative refractive index is related to the diffusive pole in the 
current correlator and that the coefficient B in (13) is related to n 2 . In figure Pt] we plot 



the ratio Re(n)/Im(n), which represents the amount of propagation with respect to the 
dissipation. Unfortunately this ratio is not large in the region of frequencies with negative 
refraction. Note that higher values of the backreaction make this ratio larger, improv- 
ing the propagation. However we do not find in the holographic superconductor a net 



improvement with respect to the case studied in [61. Moreover, as observed in 14 this 



13 The permeability (|9| is an effective permeability which is not simply related to any response function. 
However the fact that Jm[ju] < could indicate some problem in the e -fi approach. This is a delicate 
issue that goes beyond the scope of this paper. 
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numerical procedure does not allow k, 2 > .3, even if this regime is physically meaningful. 
It would be interesting to confirm our predictions for the case of larger backreaction with 
T ~ T c by an analytical computation. 

We conclude with a discussion about the validity of our results. In section [T] we 
used a long-wavelenght expansion to relate the electrical permittivity and the magnetic 
permeability to the Green function. The validity of this expansion is verified a posteriori 
by the on shell relation |k| = |n(w)|w <C 1 . We explicitly verified that this inequality 



is satisfied in every plot that we have shown. As observed in 32 if this approximation 



breaks down it is a signal that the dispersion relation has not been correctly solved. A 



detailed analysis of this case was given for the hydrodynamical diffusive pole in 33 . 

A second constraint comes from the requirement of the validity of the perturbative 
expansion in the electrical charge q. Indeed in the hydrodynamical derivation we assumed 
that the external field A ext is the dominating one. This analysis ignored the propagation 
of the induced internal field. By taking into account its contribution some constraints 
are imposed on the validity of the q 2 expansion. This requirement enforces another lower 
bound on the frequency. We checked that there are regimes of parameters in which this 
bound does not affects the whole region of negative refraction. Anyway we leave deeper 
analysis for further studies. 

A last observation is related to the fact that we consider a translation invariant system. 
We have checked that the refractive index is still negative even if we add a small amount 
of impurities that break the translation invariance. However this point requires further 
quantitative studies because it is crucially related to any concrete comparison to physical 
systems. 

Conclusions 

In this paper we have observed that the two fluid description of superconductors allows 
the existence of negative refraction in the low frequencies regime. We reached this conclu- 
sion from the analysis of the general behavior of the correlation functions of the transverse 
current correlator in the hydrodynamical regime, and we corroborated this claim by an- 
alyzing an explicit example taken from high energy physics. This example deals with 
some powerful tools developed in the context of gauge/gravity duality. Indeed the corre- 
lation functions and the transport coefficients of a general class of superconductors can 
be exactly computed from the AdS/CFT dictionary. 

On the gravity side our superconductor corresponds to an AdS/RN black hole coupled 
to a scalar field. In the broken phase this model is the dual of the EM symmetry breaking 
which leads usually to superconductivity. The EM properties of this system were already 



investigated in 16 , where it was shown that the refractive index remains positive if the 
background metric is kept fixed. Here instead we studied the full backreacted example 
and we showed that, as expected from the hydrodynamic description, negative refraction 
is allowed. 

Many interesting further lines of research are possible. A first important step should 



17 



be to set some numerical scale to the prediction of negative refraction in holographic 
superconductor, and try to make some connection with experimental physics. 

Then it would be interesting to look at the behavior of the Josephson junction: two 
layered superconductor separated by another material (an insulator, a conductor or an- 
other superconductor). They are high-T c anisotropic superconductors and are promising 
metamaterials candidates because they are easier to fabricate and less dissipative then 
usual metamaterials 34 . From the holographic perspective a gravity dual description has 



been proposed in 35 and then developed and extended in 36 38 . 

Another interesting line of investigation would be the extension of our work to the 
AdS4 case. Indeed in the normal phase it was already shown that the negative refractive 
index is allowed [7], and a similar parallelism is expected in the broken phase. 

Finally, it would be interesting to study the current correlation function for other sys- 
tems in absence of the translation invariance pole and outside the hydrodynamic regime, 
to find if they can show negative refraction. The AdS/CFT is a natural tool to study 
this problem: it can indeed describe globally charged stable systems without the first 



contribution in (13). In this case, to claim something about the refractive index of the 
system, we cannot rely on the hydrodynamical approach and we need to compute higher 
order term in the w, k expansion. 
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Figure 1: Real part of the permittivity as a function of w for k 2 
T/T c = .9 (red) .75 (blue) .6 (green) and .45 (orange). 
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Figure 2: Imaginary part of the permittivity as a function of w for k, 2 = 10 10 , .1, .2, .3 
and for T/T c = .9 (red) .75 (blue) .6 (green) and .45 (orange). 
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Figure 3: Real part of the permeability as a function of w for k 2 = 10 10 , .1, .2, .3 and for 
T/T c = .9 (red) .75 (blue) .6 (green) and .45 (orange). 
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Figure 4: Imaginary part of the permeability as a function of w for k 2 = 10 10 , .1, .2, .3 
and for T/T c = .9 (red) .75 (blue) .6 (green) and .45 (orange). 
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Figure 5: Index noi as a function of w for k, 2 
.75 (blue) .6 (green) and .45 (orange). 
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Figure 6: Zoom of the index Udl m the region of significant negative refraction (low w) 
for k 2 = KT 10 , .1, .2, .3 and for k 2 = 1(T 10 , .1, .2, .3 and for T/T c = .9 (red) .75 (blue) .6 
(green) and .45 (orange). 
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Figure 7: Ratio Re(n)/Im(n) for h 2 = KT 10 , .1, .2, .3 and for T/T c = .9 (red) .75 (blue) 
.6 (green) and .45 (orange) in the low w regime. 
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